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ABSTRACT 

We present a measurement of the if-band luminosity function (LF) of field galaxies 
obtained from near-infrared imaging of a sample of 345 galaxies selected from the 
Stromlo-APM Redshift Survey. The LF is reasonably well-ht over the ten magnitude 
range —26 < Mk < —16 by a Schechter function with parameters a = —1.16 ± 0.19, 
M* = —23.58 ± 0.42, (p* = 0.012 ± 0.008 Mpc“^, assuming a Hubble constant of 
Hq = 100 km s“^Mpc“^. We have also estimated the LF for two subsets of galaxies 
subdivided by the equivalent width of the Ha emission line at EW (Ha) = lOA. There 
is no signihcant difference in LF shape between the two samples, although there is 
a hint Icr significance) that emission line galaxies (ELGs) have M* roughly one 
magnitude fainter than non-ELGs. Contrary to the optical LF, there is no difference 
in faint-end slope a between the two samples. 

Key words: cosmology: observations — galaxies: luminosity function, mass fnnction 
— surveys 



1 INTRODUCTION 

Deep, near-infrared K-band (2.2pm) galaxy surveys are a 
powerful tool for st udyin g galaxy evolu tion, eg. Gardner, 
Cowie fc Wa inscoat (199S), Cowie et al. (1994), Glazebrook 



et al. ( 1995). Gompared to blue-optical light, near-infrared 



light is a better tracer of mass in evolved stars and the cor- 



been measured for 1797 galaxies with bj < 17.15 over a very 
large volume of space. The solid angle of the survey is 1.3 sr 
and the median redshift is about 15,300 km/s. The key to 
measuring the faint-end of the K-band LF without a huge 
investment of telescope time is to observe galaxies selected 
by their intrinsic luminosity rather than their apparent flux. 



rection for redshift dimming (the “fc-correction” ) is approx- 


cal luminosities are correlated (Mobasher, Ellis & Sharpies. 


imately 


independent of morphological type. The rapid evo- 


198C; saracco, Ghincarini & lovino 1996), in order to prefer- 



not seen in the K band, eg. Koo & Kron (1992). However, 
it is vital to have a reliable determination of the K-band 
luminosity function (LF) for nearby galaxies in order to in- 
terpret faint galaxy counts and to calculate the clustering of 
if-selected galaxy samples. Local Lf-band luminosity func- 
tions have been measured fro m op tically-selected samp les by 



Mobasher, Sharpies & Ellis (199S) and Szokoly et al. (199f ) 



and from K-select ed sam ples by Glazebrook et al. (199 



and Gardner et al. (1997). Since all of these surveys are flux- 
limited, the majority of galaxies have K-band luminosities 
close to I/)f. Even in the largest sample (Gardner et ah, 510 
galaxies), it was feasible to measure the LF over a range of 
only 5 magnitudes, and so the the faint-end slope of 

the K-band luminosity function, so important for predicting 
galaxy number counts, is still rather poorly constrained. 

In this paper we present a new estimate of the K-band 
LF over a range of ten magnitudes, based on a sub- 

sample of galaxies selected from the Stromlo-APM galaxy 
survey (Loveday et al. 1996). This survey is an ideal source 
for estimating the A-band LF since redshifts have already 



entially select galaxies of high and low luminosity and thus 
sample the luminosity range more uniformly than a flux- 
limited sample. One is thus able to measure the luminosity 
function to fainter luminosities than from a flux-limited sam- 
ple of similar size. This sampling strategy is described in M 
and the observations and data reduction are discussed in ® 
Our method of estimating ^[Mk) from a 6j-selected sam- 
ple is given in ^ and we test this method in Our results 
are presented in and we conclude in §^. For notational 
convenience, we will denote absolute magnitudes in the K 
and bj bands by Mk and Mb respectively. Throughout, we 
assume a Hubble constant of Ho = 100 km/s/Mpc. 



2 SAMPLE SELECTION 

The aim in selecting a subset of Stromlo-APM galaxies for 
which to obtain A-band photometry was to sample the mag- 
nitude range —22 < Mb < — 13 (the full range of bj abso- 
lute magnitudes in the Stromlo-APM survey) as uniformly 
as possible. An added complication in defining the sample 
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arose because we wished to obtain optical CCD images for 
the same sample of galaxies. One planned use of this op- 
tical imaging is to measure morphological parameters for 
a representative sample of galaxies at low redshift in or- 
der to compare w ith HST observations of galaxies at high 
redshift, 2 > 0.4, (Brinchmann et al. 1996). To obtain com- 
parable linear resolution to the HST data required observ- 
ing galaxies at 2 < 0.04, assuming ground-based seeing of 

1.3 arcsecond. Our primary sample thus consists of galaxies 
at redshifts 2 < 0.04. We divided the magnitude interval 
— 22 < Mb < —13 into 90 bins each of width 0.1 mag. We 
then randomly selected up to six galaxies from the Stromlo- 
APM survey with 2 < 0.04 in each bin. Due to the redshift 
limit of 2 < 0.04, this primary sample contains rather few 
galaxies brighter than Mb = —20. We therefore formed a 
supplementary sample, consisting of galaxies at 2 > 0.04 
to top up each magnitude bin, where possible, to six galax- 
ies. This supplementary sample consists entirely of galaxies 
with Mb < —20. The primary sample contains 283 galaxies, 
and the supplementary sample contains 80 galaxies, giving 
a total sample size of 363 galaxies. 



3.2 Mosaicing 

At this stage of the reduction we had ten frames per galaxy, 
with the galaxy at a different position on each frame. Mo- 
saicing of the frames and subtraction of sky background 
was performed using the dimsum p ackag e. As described by 
Stanford, Eisenhardt & Dickinson (1995), dimsum employs 



a two-pass procedure in order to mask out faint as well as 
bright images when constructing the sky background, re- 
sulting in a much flatter sky than is obtained from median 
filtering of the galaxy frames with outlier rejection. Indi- 
vidual galaxy frames were block replicated by a factor of 4 
in each dimension, allowing alignment to be performed by 
integer offsets without interpolation. Co-added images and 
corresponding exposure weight maps were made by summing 
the aligned images and were finally block averaged 2 x 2 in 
order to economize on disc space. 

For the 11 “biggrid” observations, the four offset sky 
frames were median Hltered with outlier rejection in order 
to estimate the sky background. The sky background was 
subtracted from each of the four on-source frames, which 
were then aligned and coadded. 



3 OBSERVATIONS AND DATA REDUCTION 

Imaging of the above sample of galaxies was carried out at 
the Cerro Tololo Interamerican Observatory (CTIO) 1.5m 
telescope using the CIRIM infrared array through the stan- 
dard K filter over the nine nights 1996 August 31 - Septem- 
ber 4 and 1997 October 19-22. The pixel size at //7.5 
is 1.16”, allowing most galaxies to be observed at 9 non- 
overlapping positions on the 256 x 256 array. Two frames 
were taken with the galaxy at the central position, thus 
yielding 10 frames per galaxy. Total integration time for each 
galaxy was 300 seconds. For 11 galaxies with angular size 
more than 100 arcsec, we obtained four on-source and four 
offset-sky integrations of 75 seconds each. In the following, 
these will be referred to as “biggrid” observatio ns. St andard 
stars were observed from the list of Elias et al. (1982). Dark 



frames and dome flats were taken at the start of each night. 

The infrafred frames were reduced using iraf, mosaiced 
with the dimsumQ package and image d etection and pho- 



tome try was performed using SExtractor ( Bertin fc Arnouti 

I^). 



3.1 Basic reduction 

The basic reduction process consisted of the following steps: 

(i) Non-linearity correction using the irlincor task with 
parameters cl = 0.9997, c2 = 0.0257, c3 = 0.0158 (Mike 
Keane, private communication), 

(ii) Subtraction of dark frame, 

(iii) Flathelding by dome flat, 

(iv) Masking of bad pixels. 



* DIMSUM is the Deep Infrared Mosaicing Software pack- 
age developed by Peter Eisenhardt, Mark Dickinson, Adam 
Stanford, a.nd .lohn Ward, and is available via, ftp from 
ftp: / /iraf. noao.edu/iraf/contrib/dimsumV2/dimsum. tar. Z 



3.3 Image detection and photometry 

We used SExtractor 2.0.15 to detect and measure images in 
the mosaiced frames. For both standard stars and galaxies 
we used the mag_best esti mate of ma gnitude. This yields 
a pseudo-total magnitude (Kron 198C) except in crowded 



helds, when a corrected isophotal magnitude is measured 
instead. Magnitude errors were estimated by combining in 
quadrature SExtractor’s estimate of the error from photon 
statistics and the difference between magnitudes measured 
using local and global estimates of the sky background. Of 
the selected sample of 363 galaxies, 351 were observed under 
photometric conditions and 345 yielded a A-band magni- 
tude with an estimated error of less than 0.1 mag (rms mag 
error = 0.04 mag). 



3.4 Calibration 

Since we observed only single band infrared imaging, we used 
the following simple relation to conve rt observed magnitudes 
k to standard CIT (Elias et al. 1982) magnitudes K: 



K — k -\- k(j kx^^ 



( 1 ) 



where ko is a zero-point offset, kx is the extinction coeffi- 
cient in the A-band and X is the airmass of the observation. 
We made a total of 74 standard star observations during the 
runs (an average of 8 per night) and we initially htted the pa- 
rameters ko and kx for all standard star observations com- 
bined, obtaining ko = 7.588 ± 0.020 and kx = 0.081 ± 0.018 
with rms residual magnitude error of 0.018. Six standard 
star observations with large residuals were omitted from the 
htting procedure, three of these were from the night of 1996 
August 31, which was partially non-photometric. Holding 
the zero-point term fixed at ko = 7.588 we then fitted the 
extinction coefficient separately for each night, with results 
shown in Table Galaxy magnitudes were converted to the 
CIT system using (|^ with ko = 7.588 and kx coefficients 
from the Table. Galaxies observed during non-photometric 
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Table 1. Standard star calibration data. 



Night 


Ntd 


^lel 




rms 


1996 Aug 31 


11 


3 


0.105 ±0.006 


0.017 


1996 Sep 01 


4 


1 


0.101 ±0.002 


0.003 


1996 Sep 02 


9 


0 


0.081 ± 0.004 


0.014 


1996 Sep 03 


9 


0 


0.082 ± 0.004 


0.015 


1996 Sep 04 


7 


0 


0.092 ± 0.004 


0.011 


1997 Oct 19 


9 


0 


0.083 ± 0.003 


0.010 


1997 Oct 20 


7 


1 


0.067 ± 0.004 


0.010 


1997 Oct 21 


11 


1 


0.075 ± 0.002 


0.007 


1997 Oct 22 


7 


0 


0.066 ± 0.005 


0.013 



^ Number of standard stars observed. 

^ Number of outliers deleted. 

^ Extinction coefficient in (^. 

conditions on the nights of 1996 Angust 31 and 1997 Oc- 
tober 20 were rejected: these cases were obvious from the 
rapidly varying sky background. 

3.5 Photometric repeatability 

Thirty-one galaxies, mostly of low A-band surface bright- 
ness, were observed on more than one occasion, which allows 
us to assess the repeatability of our photometry. In these 
cases, we ran SExtractor on the mosaiced galaxy frames both 
before and after coadding observations. We obtain final pho- 
tometry from the coadded frame and use the two individual 
frames and the coadded frame to estimate the rms error. 
In four cases, one observation was made using “biggrid”: in 
these cases we do not coadd the observations and instead 
use the DiMSUM-reduced observation for the final magnitude 
and estimated error. 

In Figure we plot both our estimated errors and the 
rms errors between repeated observations as a function of 
K magnitude and mean surface brightness fiK- The sym- 
bols indicate the rms magnitude error and the length of the 
line shows the estimated magnitude error from the coadded 
image. Thus if our estimated errors are a good estimate of 
the true rms then the lower ends of the error bars should 
reach zero. In most cases, the rms error is within 2-3 times 
the estimated error and in a few cases the error is slightly 
overestimated. The five points with rms > 0.4 all correspond 
to images which were broken up by SExtractor’s deblender 
(open symbols). We attempted to sum the flux from the 
components, but clearly the K magnitude for these objects 
is accurate to only ~ 0.5 mag. We see that most of the 
discrepant points are of extremely low surface brightness, 
Uk ^ 21.5 mag arcsec“^. For unbroken images, all rms er- 
rors are less than 0.2 mag and the estimated errors provide 
a reasonable estimate of the rms. The stars denote objects 
in which one of the observations used separate sky expo- 
sures (“biggrid” observations). For two out of four of these 
objects, we see an rms magnitude error larger than 0.1 but 
with a negligible estimated error (from the non-biggrid ob- 
servation). The reason for this is the poorer sky subtraction 
of the biggrid observations compared with the dimsum re- 
duction of the majority of our observations. 




8 10 12 14 16 



K 




18 20 22 
Mk 



Figure 1. Repeatability of galaxy photometry. The top panel 
plots the rms K magnitude between repeated observations against 
K magnitude from the coadded data. The symbols indicate the 
rms magnitude error and the length of the line shows the esti- 
mated magnitude error. The lower panel plots the same informa- 
tion but as a function of mean A-band surface brightness 
within the measurement aperture. Open symbols denote sources 
deblended by SExtractor, stars indicate a “biggrid” observation. 

3.6 Matching to APM 

Since the mosaiced frames cover an area of sky around 8 
arcmin on a side, albeit not to uniform depth, they pro- 
vide A-band photometry for many objects in each field in 
addition to the target galaxy. We therefore matched the im- 
ages detected by SExtractor with images in the APM scans. 
Using the matched objects in each frame, we calculated a 6- 
parameter transform from CIRIM pixel coordinates to APM 
plate coordinates and thence to RA & Dec. 

The A-band photometry for the target galaxies is pre- 
sented in the Appendix to this paper. 

3.7 Sample properties 

Since our source survey is limited by bj flux, our A-band 
sample is not a complete one, as illustrated in Figure ^ 
At faint A magnitudes, only blue objects will be in the 
Stromlo-APM sample. We are roughly complete to A « 12, 
the bluest galaxies can be seen as faint as A « 15. The 
luminosity function estimator described in corrects for 
this A-band incompleteness and allows us to use all galaxies 
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B 10 12 14 16 



K 




My. 

Figure 3. Rest-frame Mb — Mk versus Mk colour-magnitude 
plot. Plus signs represent early-type galaxies, asterisks late-type 
galaxies and dots represent unclassified galaxies. The solid line 
shows a least-squares fit to all galaxies, the dashed line a fit to 
early types and the dot-dashed line a fit to late type galaxies. 

bj magnitudes for galaxies bj ^ 15 and a few instances of 
the APM scans breaking up a very large, bright galaxy into 
fragments. An improved colour-magnitude relation will be 
available when we incorporate optical CCD magnitudes for 
these galaxies (Loveday & Lilly, in preparation). 



Figure 2. Illustration of observational selection effects in our 
sample. The upper panel plots observed bj — K colour versus 
apparent K magnitude. The line shows the apparent magnitude 
limit for the Stromlo-APM survey, bj = 17.15. The lower panel 
plots redshift against apparent K magnitude. The line in this case 
separates the primary and supplementary samples at ^ = 0.04. 



with A-band photometry. In the lower panel of this Figure 
we plot the redshift distribution as a function of apparent K 
magnitude. The requirement that z < 0.04 for the primary 
sample has a noticeable effect on the overall redshift distri- 
bution: this effect is also accounted for by our luminosity 
function estimator. 

To calculate absolute magnitudes we use galactocentric 
recession velocities and assume k-corrections in the bj band 



as given in Table 2 of Efstathiou, Ellis & Peterson (1988). 



We assume a A-band k-correction of —2.52 for all galaxy 
types (ie. the near-infrared k-correction brightens galaxies 
with redshift). This is a very good appro ximat ion to the A- 
band k-correction of Glazebrook et al. (1995) for redshifts 
2 < 0.15. 

Figure ^ shows the rest-frame {Mb — Mk) versus Mk 
colour-magnitude relation for our data. The fit to all galaxies 
is given by 



(Mb - Mk) = -0.212 x Mk - 1.20, a = 0.77. 



(2) 



The slopes in this relation for early and late-type galaxies 
(—0.135 ±0.036 and —0.209 ±0.024 respectively) are consis- 
tent with those found for E-SO (—0.095 ±0.013) and S a-Sdm 
(— 0.24±0.03) galaxies by Mobasher, Ellis & Sharpies (1986), 
although we find a significantly larger scatter about the re- 
lation. At least part of this scatter is due to saturated APM 



4 ESTIMATING THE A-BAND LUMINOSITY 
FUNCTION 



When one has a sample selected on optical bj magnitude 
and wishes to estimate the A-band luminosity function, 
the best way to proceed is to calculate a bivariate lumi- 
nosity function (BLF, cf){MK , Mb)) allowing for known se- 
lection in bj flux and Mb absolute magnitude and then 
to integrate over Mb to obtain 4>{Mk)- One can estimate 
the shape of 4>{Mk, Mb), independently of inhomogeneities 
in the galaxy distribution, using the m aximu m likelihood 
method of Sandage, Tammann & Yahil ( 197S| ). The proba- 
bility of seeing a galaxy with A-band luminosity Mk and B 
band luminosity at redshift Zi in our sample is given by 



Pi 



<!>{Mk,L],)S{L],) 



L 






Mi) 



4>{LK,LB)S{LB)dLKdLB 



( 3 ) 



The function S{Lb) accounts for the known selection in ab- 
solute B magnitude and the luminosity limits Lsminizi) and 
LBuiaxi^i) are the minimum and maximum A-band lumi- 
nosities observable at redshift Zi in a sample limited by ap- 
parent b magnitude. For the sample analysed here, there are 
no flux limits in the fc-band, and so the integral over A-band 
luminosity runs from 0 to ±oo. The maximum- likelihood 
shape of the BLF (j){MK , Mb) is estimated by maximizing 
the likelihood C — Pi (the product of the individual 

probabilities Pi for the Ng galaxies in the sample) with re- 
spect to the parameters describing the BLF. 

In practice, we do not have a good a priori para- 
metric model for ({/{Mk , Mb), and so instead we measure 
<j){MK, Mb) in a non-parametric way using an extension 
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Figure 4. Sample selection as a function of Mb absolute magni- 
tude for the primary (thick lines) and supplementary (thin lines) 
samples. The upper panel shows histograms of Mb for all 1797 
galaxies in the Stromlo-APM survey. The middle panel shows the 
same for the 345 selected galaxies with Tf-band photometry. The 
lower panel shows the sampling functions S{Mb), the ratio of 
selected to total galaxies in each bin. 



of the Efstathiou, Ellis & Peterson 



stepwise maxi- 



mum likelihood (SWML) method. Sodre & Lahav (199f) 
have extended the SWML method to estimate the bivariate 
diameter-lnminosity function and to allow for sample incom- 
pleteness. We adopt their extension of the SWML estima- 
tor here, inclnding the sampling fnnction S{Mb) separately 
for the primary and supplementary galaxy samples. These 
sampling fnnctions are illustrated in Figure which shows 
histograms oi Mb for the entire Stromlo-APM sample, for 
the selected galaxies with good if-band photometry and the 
sampling functions S{Mb)- 

We measure the bivariate luminosity function in bins of 
absolute K and B magnitude. 



4>{Mk,Mb) = (jijk, j = 1, . . . , fe = 1, . . . , (4) 

where Mf — AMkI‘2 < Mk < Mf -I- AMk/‘2 and Mb — 
AMb/2 <Mb <M%+ AMb/2. 

The log-likelihood is given by 

Ng 

ln£ = ^ y] y] Wi,fcln[,)),fe5(M|)l- 

i=l j = \ k = l 



^^9 

El- 



EE 

3 = 1 k = l 



-|- const. 



(5) 



Here Wijk = W(Mf - Mf, M%- Ml), with, 

r 1 if - AMk /2 < a; < AMx /2 
W{x, y) = and - AMs/2 <y< AMb/2, (6) 

[ 0 otherwise. 

The sampling fnnction S{Mb) (Fig. ^) is incorporated into 
the ramp function H . Writing = Mf — AMk/ 2, Mf = 

Mf + AMk/ 2, M~ = Ml- AMs/2, and M+ = Ml + 
AMb /2, then 



cM' 



I'M' 



Bijk — 



AMk AMb 



UMk / S{MB)dMB, 



(7) 



where Mf = max[M^, min(M)^, MKiim)] a^nd Ml = 
max[M^, min(M^, Msum)]. 

To fix the otherwise arbitrary normalisation constant 
we apply the constraint 



<7 = E E 

3 k 



AMk AMb -1 = 0 

Lf Lf J 



(8) 



with (5 = 1.5 and the fiducial luminosity Lf corresponding 
to M — —20 using a Lagrangian multiplier A. The new like- 
lihood \nC' = \nC + Xg is maximised with respect to the 
(fijk and A, requiring that A = 0. 

Setting din /dipjk = 0, one arrives at a maximum 
likelihood estimate for 4>jk, 

HifkAMKAMs 
^ Em=f <t>lmHumAMKAMB') 

where rijk = X/ilfi 4he number of galaxies in the 

(/, k)th bin. Errors in tfijk are estimated via the inverse of 
the information matrix (Efstathiou et al. 1988). 

As with all density-independent estimators, informa- 
tion about the overall normalisation is lost. We therefore 
normalise our BLF to the mean density of galaxies with 
—22 < Mb < —13 in the full Stromlo-APM sample, h = 
0.071h ^ Mpc~^, calculated as described by Loveday et al. 
( 1992|) . (Note that the density h = 0.047/i®Mpc“® quoted by 
Loveday et al. (1992) is for the restricted magnitude range 
-22 < Mb < -15.) 

Finally, we obtain the Tf-band luminosity function by 
summing over B magnitude bins, 

mi) = E 

fc=i 



jk — k 







One c an then fit a given functional form, eg. a Schechter 
(1976|) function, to the stepwise estimate by least-squares. 



5 TEST OF THE METHOD 



We have tested the above procedure by using it to estimate 
the A'-band luminosity function from a set of Monte Carlo 
simulations. We ge nerate d nine mock Stromlo surveys by a 
Soneira & Peebles (197f ) hierarchical clustering simulation 
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Figure 5. Test of our method to estimate the it'-band luminosity 
function from fej-limited samples. The upper panels show the K- 
band LF estimated from each of nine Soneira-Peebles simulations 
and the lower panel shows the mean and rms (error bars offset to 
right). In each case the curve shows the input LF. 



with similar clustering properties to that mea sured from the 
Stromlo-APM Survey (Loveday et al. 1995). Each galaxy 
in the simulation was assigned a if-band luminosity drawn 
at random from a Schechter function with a = —1.20 and 

= —23.6 and then assigned a bj magnitude according 
to our observed colour-luminosity relation (^) . Galaxies were 
selected on their apparent bj magnitude, bj < 17.15. This 
process was repeated until each simulation contained 2000 
galaxies. We then sampled each simulation by absolute Mb 
magnitude as described in §^, Hnally yielding an average of 
359 galaxies per simulation. We calculated the K-band lumi- 
nosity function 4>{Mk) for each simulation as described in 
and fit a Schechter function to each by least squares. Aver- 
aging over the nine simulations, and estimating the BLF in 
bins of width 0.5 mag, we measure mean and rms Schechter 
function parameters a = —1.17 ± 0.07, M* = —23.6 ± 0.2. 
The errors on the mean values are times smaller than 
the quoted rms scatter between the simulations. Our esti- 
mates of a and M* are in excellent agreement with the input 
luminosity function. 

The SWML estimates of the Js'-band LF for the simu- 
lations are shown in Figure The points show the SWML 
estimate of 4>{Mk) from the nine simulations and the curve 
shows the input Schechter function with shape a = —1.20, 

= —23.6. The error bars going through the data points 
show the errors determined from the covariance matrix. The 



lower panel shows the average over the 9 simulations. The 
rms scatter between realisations for each data point (shown 
by the error bars offset slightly to the right) are in satisfac- 
tory agreement with the predicted errors. 

Overall, we find that our procedure for estimating 
4>{Mk) from a sample limited by apparent b magnitude and 
further selected by absolute B magnitude provides a robust 
and unbiased estimate of the Lf-band LF over a wide range 
of absolute magnitudes. 



6 RESULTS 



Our estimated TL-band luminosity function is shown in Fig- 
ure^ The curve shows a Schechter function ht to the SWML 
estimate using least squares. We allow for Hnite bin width 
by calculating the mean predicted {(fij) at the absolute mag- 
nitude of each galaxy in each bin, rather than simply cal- 
culating (fi{M) at the bin centre. The inset shows 1 and 
2 a likelihood contours for the shape parameters a and 
M^, where the normalisation (j>* is adjusted to maximize 
the likelihood at each grid point. The best-ht Schechter pa- 
rameters are a = —1.16 ± 0.19, = —23.58 ± 0.42 and 

(\>* — 0.0121 ± 0.0082/i^Mpc“®. The quoted errors on a and 

come from the bounding-box of the Icr likelihood con- 
tour. 

Clearly the characteristic magnitude MJ is rather 
poorly constrained from these data, and so in order to mea- 
sure the faint-end slope a independent of M^, we have also 
ht a straight line to our estimated LF over the restricted lu- 
minosity range M > —22. We measure a slope — 1.08T0.15, 
consistent with that measured from the Schechter ht. Fur- 
thermore, this slope does not change signihcantly varying 
the bright luminosity cut from —20.5 to —23. 

Since we normalised the LF to the same number density 
as the full Stromlo-APM sample, the error in (f)* is domi- 
nated by the uncertainty in shape of the LF. The Schechter 
function provides a reasonable ht to the SWML estimate 
over the full range of ten magnitudes; there is no indication 
of any faint-end turnup. 

In Figure ^ and Table |^, we compare our new estimate 
of the TL-band luminosity function with those of other work- 
ers. Following Glazebrook et al. (1995) and Gardner et al. 
(1997), we have added 0.22 mag to the points from Figure 2 
of Mobasher et al. (1993) to account for their method of 
calculating K-corrections. We have also applied an aperture 
correction of —0.30 mag to the Glazebrook et al. data (see 
Gardner et al. 1997). Ours is the hrst estimate of 0(Mx) 
fainter than Mk = —20. Brighter than this all estimates are 
in reasonable agreement, particularly given the uncertain- 
ties in the normalisation of the LF from small samples of 
galaxies. The normalisation of the Glazebrook et al. (1995) 
sample is about twice that of all the other samples. This is 
likely to be due to sampling fluctuations, since their survey 
covers only 552 arcmin^ and contains a total of only 124 
galaxies. 

Note that the bright-end points of our estimated LF 
have drop ped s ignihcantly since the preliminary analysis of 
Loveday (1998). This is largely due to the improved K-band 
photometry obtained by mosaiccing the galaxy frames with 
DIMSUM compared with a cruder algorithm used earlier. 
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Table 2. Schechter function fits to it'-band LFs. 



Sample 


a. 




M 


* 

K 




4>* jh' 


-3 


Mpc^ 


Mobasher et al. 1993 


-1.0 


± 


0.3 


-23.4 


± 


0.3 


0.0112 


± 


0.0016 


Glazebrook et al. 1995 {z < 0.2) 


-1.04 


± 


0.31 


-23.02 


± 


0.23 


0.029 


± 


0.007 


Gardner et al. 1997 


-0.91 


± 


0.24 


-23.12 


± 


0.17 


0.1 


0166 


Szokoly et al. 1998 


-1.3 


± 


0.2 


-23.6 


± 


0.3 


0.012 


± 


0.004 


This work (all galaxies) 


-1.16 


± 


0.19 


-23.58 


± 


0.42 


0.012 


± 


0.008 


This work (non-ELG) 


-1.3 


± 


0.4 


-23.5 


± 


0.7 


0.001 


± 


0.001 


This work (ELG) 


-1.2 


± 


0.4 


-22.5 


± 


0.8 


0.011 


± 


0.013 





Figure 8. Comparison of the JC-band luminosity function esti- 
mated from our non-ELG sample (filled symbols, continuous line) 
and ELG sample (open symbols, dashed line). 



6.1 ELG versus non-ELG 

As an aid to understanding the processes of galaxy forma- 
tion and evolution, it is of great help to have estimates of the 
luminosity function for different types of galaxies. The most 
objective way of separating Stromlo-APM g alaxies is via the 
equivalent width of the Ha emission line (Loveday, Tresse 
& Maddox 199£, hereafter LTM). Due to the small size of 
our if-band sample, we subdivide into just two subsamples, 
those with EW (Ha) < lOA (non-ELG, 138 galaxies) and 
those with EW (Ha) > lOA (ELG, 134 galaxies). For the re- 
maining 71 galaxies, no EW (Ha) measurement is available 
(LTM). For both subsamples we recalculated the sampling 
function S{Mb) appropriate for galaxies with EW (Ha) less 
or greater than lOA as appropriate. The LFs were nor- 
malised to mean densities of h = 0.018 h^Mpc”^ (non-ELG) 



and h = 0.098 h^Mpc”^ (ELG) as determined for galaxies 
with these emission line properties from the full Stromlo- 
APM sample. The fact that the estimated mean density 
of ELGs is larger than that for galaxies of any EW (Ha) 
{n = 0.071h®Mpc“^, cf. ^ is a consequence of the mean 
density of subsamples of Stromlo-APM galaxies being de- 
termined to only ~ 30% accuracy (LTM). 

The estimated LFs for each subsample are shown in Fig- 
ure^ and the la likelihood contours for the best fit Schechter 
parameters are shown in Figure The Schechter function 
parameters themselves are listed in Table Compared with 
the optical luminosity functions for EW (Ha)-selected sam- 
ples (LTM), the Js'-band LFs for ELGs and non-ELGs are 
very similar in shape, with just a ~ Icr indication that ELGs 



Figure 9. la likelihood contours for the best fit Schechter pa- 
rameters for the non-ELG and ELG samples. 



have characteristic magnitudes M* about one magnitude 
fainter than non-ELGs. Note that the faint-end slopes for 
the two samples are completely consistent. These results are 
in accord with the LFs estimated for E/SO and spiral galax- 
ies by Mobasher et al. (1993). 



7 CONCLUSIONS 

We have presented a new estimate of the A'-band luminosity 
function obtained from infrared array imaging of a subsam- 
ple of galaxies from the Stromlo-APM survey. We measure 
4>{Mk) over a range of 10 magnitudes, a significantly 

greater range of luminosities than has been measured be- 
fore now. Our LF is consistent with earlier estimates, and 
shows that a Schechter function of moderate faint-end slope 
(a = —1.16) provides a good fit to the LF as faint as M* -1-8 
with no hint of any upturn at the faint-end. This observa- 
tion appears to rule out the suggestion of Glazebrook et al. 
(1995) that there may be an extra dwarf component to the 
K-hand LF, or at least that if there is such a component, 
it must have characteristic magnitude Mdwarf ^ - f 8. In 
contrast, the faint-end upturn seen in the field bj LF (Love- 



day 1997) is apparent by Mg -|- 6, and in the CfA Redshift 
Survey ( [Marzke, Huchra fc Geller 1994 ) by -|-3. In the 
Coma cluster, an upturn in the H-hand (1.6u.m) LF is also 
seen at Mg + 3 ( De Propris et al. 1998 ). Note that we do 
see a slight upturn in (^{Mk) near Mg + 3, faintward of the 
three low points around Mg = —22, but the LF flattens out 
again by Mg = — 20. 

The TF-band LFs for non-ELGs and ELGs are very sim- 
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ilar in shape, with no difference in faint-end slopes, contrary 
to what is seen in the optical LF (LTM). The steep optical 
LF for ELGs might be co nsistent with a ~ fiat ma ss function 
in a fading burst model (Hogg & Phinney 1997), although 
the particular model that Hogg & Phinney consider would 
steepen the K-hand LF even more than the optical LF. 

The main limitation with our sample is that the bj se- 
lection of the source survey introduces a bias against red 
galaxies compared to a TL-selected sample. However, the 
good agreement between the TL-band LFs determined from 
K and fej-limited samples in Fig. ^ suggests that our bj se- 
lection does not significantly affect the estimated K-hand 
LF. In particular, it is unlikely that we are underestimat- 
ing the faint-end slope of the Ji'-band LF since red galax- 
ies tend to be luminous in K (Mobasher^ Ellis fc Sharpies 



198C; Saracco, Chincarini & lovino 1996). We will further 
address the issue of colour selection in a future paper in 
which we plan to parameterize the bivariate LF 4>{Mk, Mb) 
using CCD bj magnitudes in addition to the 7L-band data 
presented here. 
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APPENDIX A: THE K-BAND DATA 



In Table Al we present a sampling of our K-band photo- 
metric data. The complete catal ogue will be available from 
the Astronomical Data Centre (http://adc.gsfc.nasa.gov/). 
The first five columns of this table come from the Stromlo- 
APM survey (Loveday et al. 1996). The subsequent seven 
columns are derived from the K-band data. Each column in 
the table is described below. 

(1) Name: Galaxy naming follows the same conv ention as 
the APM Bright Galaxy Catalogue ( Loveday 1996| ) and the 
Stromlo-APM Redshift Survey (Loveday et al. 1996). 

(2) , (3) RA, dec: Right ascension (hours, minutes, sec- 
onds) and declination (degrees, arcminutes, arcseconds) in 
1950 coordinates. 

(4) 6j: bj magnitude. 

(5) cz'. Heliocentric recession velocity in km/s. 

(6) , (7) K, Kerr: K magnitude and its estimated error. 

(8), (9) Maj, Min: Semi-major and minor axes of measure- 
ment ellipse in arcseconds. 

(10) PA: Position angle in degrees measured clockwise 
from south-north line. 



(11) Flags: Flags output by SExtractor (Bertin fc Arnout; 



1996|). 



(12) Notes: 1: Galaxy deblended by SExtractor; flux of 

components was summed. 

2: “Biggrid” observation: photometry may be less reli- 

able than most galaxies. 



NGC1672: APM deblended this galaxy. 
Possibly non-photometric. 

Central galaxy in close group of 3. 
Double nucleus. 
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Table Al. / C-band galaxy photometry . 



Name 




RA 




Dec 


bj 


cz 


K 


^^err 


Maj 


Min 


PA 


Flags 


Notes 


075 + 069-077 


21 


26 


16.37 


-68 


39 


53.3 


16.32 


11033 


12.50 


0.02 


16 


11 


46 


3 




076 - 113-015 


22 


56 


20.64 


-69 


43 


06.6 


16.48 


3813 


13.25 


0.04 


32 


14 


61 


2 




077 + 055+032 


23 


11 


33.04 


-70 


38 


10.8 


17.05 


33310 


14.09 


0.05 


13 


8 


104 


3 




077 + 062-116 


23 


11 


36.12 


-67 


52 


12.9 


16.57 


29923 


11.98 


0.02 


19 


15 


172 


3 




078 - 130-118 


00 


21 


06.17 


-62 


46 


46.8 


16.52 


12128 


13.55 


0.04 


12 


7 


116 


3 


1 


078 + 109+066 


23 


39 


46.85 


-66 


14 


01.5 


14.62 


10150 


11.91 


0.07 


44 


20 


86 


2 


2 


078 + 012+091 


23 


57 


37.39 


-66 


47 


31.6 


15.42 


21754 


11.05 


0.01 


25 


23 


69 


18 




080 - 018-033 


01 


31 


33.11 


-64 


28 


28.6 


17.01 


8088 


13.95 


0.05 


16 


12 


136 


0 




080 + 009+024 


01 


26 


43.01 


-65 


31 


49.0 


14.84 


1624 


11.47 


0.01 


25 


22 


12 


0 




082 + 032-078 


02 


50 


37.37 


-63 


38 


30.4 


16.68 


30356 


13.16 


0.03 


15 


7 


110 


2 




107 - 092-049 


21 


31 


44.58 


-64 


06 


47.1 


16.14 


3179 


8.54 


0.00 


86 


43 


174 


48 




107 - 009+013 


21 


17 


34.35 


-65 


19 


12.8 


16.66 


5118 


14.06 


0.07 


18 


16 


175 


2 




107 - 053-114 


21 


24 


40.36 


-62 


56 


09.3 


17.09 


8507 


14.49 


0.08 


19 


8 


152 


0 




108 + 093+006 


21 


43 


24.76 


-65 


10 


25.7 


16.03 


10497 


13.20 


0.04 


18 


14 


24 


0 




108 - 114+028 


22 


20 


24.21 


-65 


32 


46.0 


16.63 


6146 


13.69 


0.05 


26 


14 


89 


0 




108 - 105+055 


22 


19 


19.57 


-66 


03 


24.1 


14.56 


10779 


11.03 


0.01 


47 


21 


47 


16 




108 - 083-130 


22 


13 


20.13 


-62 


38 


40.3 


16.85 


36934 


12.85 


0.03 


14 


13 


123 


0 




109 + 063+033 


22 


32 


40.84 


-65 


41 


37.3 


16.69 


21798 


11.97 


0.02 


18 


14 


130 


0 




109 + 014+028 


22 


41 


34.01 


-65 


37 


16.2 


14.84 


3269 


11.73 


0.02 


26 


25 


31 


0 




110 + 020+024 


23 


24 


18.95 


-65 


32 


45.2 


15.04 


1991 


11.30 


0.02 


82 


38 


150 


18 




111 + 063-007 


23 


50 


48.93 


-59 


58 


26.5 


14.81 


3320 


13.05 


0.54 


53 


25 


31 


0 


2 


111 + 118+005 


23 


42 


33.53 


-60 


08 


39.0 


15.52 


3414 


14.24 


0.24 


30 


18 


73 


2 


2 


111 + 025+058 


23 


56 


29.36 


-61 


11 


27.6 


16.20 


28872 


12.24 


0.02 


24 


15 


179 


0 




112 + 092-034 


00 


24 


46.45 


-59 


24 


30.9 


15.97 


11668 


12.15 


0.02 


18 


12 


75 


0 




112 - 083-068 


00 


50 


12.34 


-58 


47 


36.5 


16.15 


5150 


14.75 


0.08 


13 


9 


180 


0 




112 + 079+009 


00 


26 


25.33 


-60 


13 


13.5 


16.62 


4720 


11.89 


0.02 


38 


27 


3 


16 




113 + 004+088 


01 


15 


15.04 


-61 


43 


42.3 


15.31 


8591 


11.54 


0.01 


30 


11 


7 


16 




114 - 031+009 


01 


58 


39.99 


-60 


14 


34.1 


17.06 


6743 


13.53 


0.04 


16 


6 


123 


0 




114 - 015+074 


01 


56 


18.04 


-61 


27 


21.5 


14.80 


7002 


9.98 


0.01 


49 


39 


108 


16 




114 - 123-003 


02 


12 


22.81 


-59 


56 


14.5 


16.35 


1471 


15.60 


0.14 


9 


5 


93 


0 




116 - 048-060 


03 


17 


06.64 


-58 


58 


03.9 


16.62 


21466 


12.84 


0.03 


11 


11 


68 


0 




116 - 074+000 


03 


21 


19.26 


-60 


04 


51.6 


17.04 


5664 


15.56 


0.09 


16 


11 


85 


3 




117 + 080-131 


03 


36 


49.50 


-57 


36 


19.7 


16.02 


4952 


12.22 


0.02 


23 


9 


31 


0 




117 + 097-087 


03 


34 


07.28 


-58 


24 


46.6 


16.11 


17774 


12.49 


0.02 


16 


8 


96 


0 




118 - 094-066 


04 


39 


27.28 


-58 


50 


19.6 


14.55 


1176 


10.64 


0.01 


26 


25 


36 


16 




118 - 129-038 


04 


44 


42.52 


-59 


18 


46.7 


15.19 


1278 


7.49 


0.00 


108 


54 


90 


16 


3 


118 - 080+121 


04 


38 


41.93 


-62 


20 


35.2 


16.76 


8488 


14.17 


0.06 


16 


9 


72 


3 


1 


118 - 094+123 


04 


41 


03.91 


-62 


21 


21.0 


16.00 


6157 


13.90 


0.05 


21 


18 


60 


0 




144 + 044-118 


20 


47 


21.43 


-57 


52 


35.6 


15.32 


3233 


12.43 


0.02 


21 


17 


155 


2 




144 - 126-031 


21 


12 


03.18 


-59 


26 


26.9 


16.24 


9493 


11.89 


0.02 


11 


11 


144 


0 




144 + 095-133 


20 


40 


26.42 


-57 


34 


32.3 


16.93 


10964 


12.52 


0.02 


16 


10 


128 


0 




144 + 091+082 


20 


39 


21.60 


-61 


34 


36.5 


15.69 


22306 


11.88 


0.02 


22 


15 


88 


2 




145 + 055+050 


21 


23 


44.04 


-61 


02 


29.3 


15.32 


4404 


10.30 


0.01 


34 


12 


82 


2 




145 - 099-026 


21 


46 


53.28 


-59 


35 


26.3 


16.36 


8053 


12.11 


0.02 


19 


14 


40 


2 




145 + 035+005 


21 


26 


58.51 


-60 


13 


18.9 


15.07 


8660 


10.61 


0.01 


43 


32 


136 


18 




147 - 070+099 


22 


59 


03.11 


-61 


53 


25.3 


17.02 


7787 


13.94 


0.05 


23 


20 


85 


2 




148 + 099-068 


23 


11 


40.26 


-58 


46 


11.3 


16.55 


3376 


13.93 


0.05 


16 


8 


30 


0 




149 - 013-101 


00 


01 


43.21 


-53 


11 


47.1 


14.60 


9773 


10.13 


0.01 


32 


25 


169 


16 




149 + 031+040 


23 


55 


51.03 


-55 


48 


43.0 


16.27 


9477 


11.62 


0.01 


22 


7 


9 


0 
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